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ABSTRACT 

The antecedent precipitation index (API) has been a useful indicator of soil moisture con- 

• • • 
ditions for watershed runoff calculations and recent attempts to correlate this index with space- 

borne microwave observations have been fairly successful. We show that the prognostic equation 
for soil moisture used in some of the atmospheric general circulation models (GCM) together 
with Thomthwaite-Mather parameterization of actual evapotranspiration leads to API equations. 
The recession coefficient for API is found to depend on climatic factors through potential evapo- 
transpiration and on soil vtexture through the field capacity and the permanent wilting point. 
Climatological data for Wisconsm together with a recently developed model for global insolation 
are used to simulate the annual trend of the recession coefficient. Good quantitative agreem.^nt 
is shown with the observed trend at Fennimore and Colby watersheds in Wisconsin. This study 
suggests that API could be a unifying vocabulary for watershed and atmospheric general circula- 
tion modelers. 
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A SIMULATION STUDY OF THE RECESSION COEFFICIENT FOR ANTECEDENT 

PRECIPITATION INDEX 

1.0 INTRODUCTION 

Soil moisture condition just before a rainfall has long been recognized as a msyor factor In 
watershed runoff predictions, although actual runoff depends upon many other factors, e.g., 
physical*, litnd features and vegetation. It is readily understood that when the soil is wet, run* 
off is greater than when it is dry. Soil moisture condition is also a major consideration in 
general circulation models (GCM) used in global models of the atmosphere because over the 
land the atmosphere exchanges heat with the surface and derives atmospheric moisture at a 
rate depending upon soil we''.ness. Local weather such as sea breeze and heat islands are also 
known to depend upon the state of the. soil moisture conditions. 

Since conventional methods for soil moisture estimation are impractical for large areas on 

a timely basis, hydrologists (see Linsley et d., 1949) ‘indexed’ the antecedent precipitation as 

a means of estimating the moisture conditions in a watershed. Tlrese index values can as well 

• . 

be applied to the large area represented by grid cells used as input locations for global models 
of the atmosphere. ■ Cells used in a GCM are generally larger than most watersheddi, making 
the estimation of the average soil moisture even more difficult when using conventional mea- 
surements. The antecedent precipitation index (API) suggested by Lindsley has been correlated 
to microwave measurements made from space platforms (McFarland, 1976 and Blanchard et d., 
1981) thu^ offering the possibility of repeat me asurements on a timely basis that could be 
readily converted to input for global models. Before extensive use of the API index, we need 
to determine the relation between soil moisture input in GOsI and the index. 

Gearly, soil left to itself will keep drying, and it is precipitation which can alter the 

. • ♦ * 

moisture conditions to any significant degree when large areas are considered. Saxton and 
Lenz (1967) expressed the API in the following recurrence form, 
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APIj - K(APIj.i Rj.i) (1) 

where K is the rece^on coefficient and Rj. | is the amount of rainfall on 0 - l)th day. Once 
the recession coefficient and an initial value of API are given, the API equation forms a pow- 

I 

erful model for simulating the moisture ^conditions. Tlie model to date remains conceptual, 
since no formal basis for this equation has been provided. Starting with the prognostic equa- 
tion for soil moisture used in some of the atmospheric GCMs we derive the API equation. 
Using climatologicai and soil texture data we siimulate the recession coefficient and show 
agreement with observations. 

2. API MODEL 

For horizontally-homogeneous vertically-stratified bare soils with known soil physical char- 
acteristics cne can use the Philip-deVries equations (Philip and deVries, 1957) to simulate the 
moisture conditions. The numerical methods for the solution of these> equations are complex, 
and becomes computationally quite expensive if ‘exact’ solutions are sought. Even accepting 
the inherent assumptions of these equations (e.g., no hystcrisis), it becomes ratlier difficult to 
provide the required basic data for the solution when areas of a few hectares and larger are 
considered. We realize point measurement and calculation of soil moisture could be exact, 
but for large areas which may contain moisture impermeable areas, depressions and vegetation 
of varied type and density the meaning of ‘exact soil moisture’ becomes uncle/iU (Engman, 
1981). It is probably more appropriate to describe the moisture conditions through some 
statistical attributes. The API equation implicitly describes Ute moisture conditions as a 
stochastic variable following the first-order markov process (see discussion below). 

In formulating the atmospheric GCM, Manabe (1969) considered the following prognostic 
equation for moisture conditions, 

Z « P(t) - Ej(t) 


( 2 ) 



where Z is the thickness of a fairly deep soil surface zone, 0 is volumetric moisture, P(t) and 

Eg(t) are, respectively, the rates of precipitation and actual evapotranspiration. Some other 

♦ » 

GCMs (e.g., Washington and Williamson, 1977) also include this soil moisture equation. It is 
essentially a one-layer moisture budget equation where the sub-layer remains passive with re- 
gard to moisture dynamics. 

There is no rigid prescription for choosing the thickness of the surface layer, Z. For 
bare soils a thickness of 0.2m may be appropriate (Sellers, 1965; Jackson et ai, 1973), but 
for vegetated soils moisture extraction may occur from a deeper soil layer. (Note tlrat actual 
root depth vary with plant species, climatic factors and soil factors.) Manabe (1969) considers 
the maximum water holding capacity of the surface layer, Z • O.lSm (df^ being the field 
capacity), which for Colby (0fg “ 0.482) and Fennimore (dfg * 0.364) watersheds to be stud- 
ied here correspond to Z values of 0.32 ant! 0.42m Noting tliat the analysis of recession co- 

* 

efficients for these watersheds is based on the soil moisture in the surface 0.31m layer, & com- 
promized value of Z * 0.35m is thouglit to be reasonable. (The implication of choosing a 
fixed value of Z rather than Z is that the rate of drying would vary with soil texture; 
sandy soils will dry faster than clayey soils. If Z 0^^ is fi.xed then tlie rate of drying for all 
soils will t? equal, but the thickness of the soil layer which is drying would depend upon soil 
tj'pe. We have opted to fix Z in order to make comparisons of predicted rate of drying with 
Saxton and Lenz’ observed rate of drying). 

Actual evapotranspiration depends upon evaporative demand of the atmosphere, transpira- 
tion through tlie plants and soil water availability. Although the Penman equation (Penman, 
1948) and some of its derivatives (e.g., Prisstley and Taylor, 1972; Jury and Tanner, 1975; 
Tliom and Oliver, 1977) are found to give fairly^ good estimates of potential evapotranspiration, 
no such equation has yet been found for actual evapotranspiration. Several linear and non- 
linear relations between soil moisture and the ratio of actual and potential evapotranspiration 
are documented (Budyko, 1956; Thomthwaite and Mather, 1957; Holmes and Robertson, 1963; 
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Eagleman, 1971; Davies and Allen, 1973; Barton, 1979 and Marsh et al„ 1981). ‘Whether any 
one of these relations is universally applicable is not known, but noting the simplicity and the 
assessment of Lowry (1959), Holmes and Robertson (1963) and Yaron et al (1973), the Thomth- 
waite-Mather relation is chosen for this study 
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where is moisture at the permanent wilting point and Ep is potential evapotranspiration. 


With eqn.(3), the solution of eqn.(2) can be written as 

t + At 

W(t + At) - ^(t + At, t) W(t) + f ^(t + At, t') P(t‘)dt' 

■'t 


where 


W(t) » Z[6(t)-e^. ] 


0(t + At, t) * exp 




E^(t')dt' 




} 


(4) 


At * time interval 

Note that if one postuates that the moisture conditions will vary stochastically according to a 
first order markov process then eqn.(4) can be written without going through the derivation. 
Our derivation of this equation considers the underlying physics and assumptions. 


If time interval At is chosen as one day and j is used as the day index then eqn.(4) can be 
expressed in several alternate forms, 

(i) If during the previous day there was no rain, 


where 


Wj = 


K W 


j- 1 


(5a) 


K “ exp (• 


E 

Z (6fc - 0^) 


) 


( 6 ) 


4 



E * daily total potential evapotranspiration for the previous day. 


(ii) If rain during previous day occurred after sunset (See Baier and Robertson, 1966) 

Wj « K Wj . 1 + Rj . I (5b) 

where Rj _ ^ is the amount of rainfall. (Note that evapotranspiration occurs mostly during the 
sunshine hours, and it is practically zero during the rain). 


(iii) If rain intermittently occurred during the previous sunshine period (compare with eqn, 

♦ 

(D). 


Wj * K(Wj.i + Rj.i) 


(5c) 


Clearly, any one of the above equations would be applicable depending upon the nature of 
the rainfall occurrance. The choice of a particular form for use throughout the year should be 
made only after analyzing the rainfall statistics of the location If rainfall statistics show that 
the occurence of rain is strongly baised toward night time hours then (Sb) may be used, and if 

there is no bais then (5c) would be more appropriate. 

• 

3. SIMULATION OF RECESSION COEFFICIENT 

The recession coefficient as given by eqn.(6) can be simulated knowing soil texture and 
climatic data. From Saxton and Lenz (1967) the maximum available water, Z(6fg - 0^), is 
0.0735m and 0.105m respectively in the top 0.35m soil layer for Fennimore and Colby water- 
sheds. Climatological data, actually observed for these watersheds are unavailable, and were 
therefore synthesized (Table 1) from published sources (Reitan, 1960; Kung et al„ 1964; Flowers 
et al, 1969; Bryson and Hare, 1974, and Climatological data of the U.S. Weather Bureau). Rec- 
ognizing the inherent data variability and the associated errors in calculating the potential evapo- 
transpiration, a consistency check will be discussed. 

• * 


Global insolation is a major factor determining potential evapotranspiration, and the quality 
of Weather Service data for insolation is questionable (Hoyt, 1978). We used a model developed 
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by Choudhuiy (1981) to calculate the daily total insolation (see Appendix). A comparative illus* 
tration of the model for daily total insolation at Rockville, MD, is shown in Fig. 1. The pertin* 
ent climatological data (atmospheric precipj table 'imter, turbidity, cloud fraction and optical 
thickness and surface albedo) used to simulate insolations at Wisconsin are given in Table 1 and • 
the calailated insolations are shown hi Fig. 2. 


The daily total potential evapotronspiration is calculated using Penman equation (Thom and 
OUver, 1977) 


E • 


L (1 t t)p 


2.6 X 10’ *^ r(e, -e) (1 + 0.54U) 

(TTT) ” 


(m) 


(7) 


where 


r- 1.192 X 10-7 (JL)t2 

- . 1 .24 TC* 

R, • (1 - a) S -*.98 X 10-3 e, T* (1 - e.) (1 - ) 


■ 6,11 exp 


17.27 (T - 273.2) 
T - 35.86 


6a - 0.7 + 5.95 X 10"^ e exp ( 


L is the latent heat of evaporation (2.47 X 10^ J kg"*), p is tlie density of water (1000 kg m"^), 
e is vapor pressure (mbar), T is air temperature (K), U is wind speed (m/sec), a is surface albedo, 
P is surface air pressure (985 mbar), R„ is net radiation 0 m"“), surface emissivity (0.97), S 
is global insolation (Appendix), e, atmospheric emissivity (Idso, 1981), c is cloud fraction and r 
is cloud optical thickness. 


As a check of the climatological data and model equations, the simulated potential evapo- 

transpiration and net radiation are shown in Fig. 3 togetlier with the regression equation based on 

♦ , 

observed evapotronspiration and net radiation during July through September at Hancock, Wiscon- 
sin (Tanner and Pelton, 1960). Although simulated values are more-or-less within the standard 
error of estimate for the regression equation, we see that the calculated evapotronspiration are gen- 
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gencirally lower than what would be expected from the regression equation. Tanner and Pelton 
(1960) also found that the Penman equation somewhat underestimates the evapotranspiraflon. 

Tlie simulated annual trend of the recession coefficient for Fennimore wateished together 
with observations (Saxton and Lenz, 1967) is shown in Fig. 4 using both Penman equation and 

t 

the regression equation from Tanner and Pelton (1960). Except for April-May period, the simu- 
lated trend is within the observed values. If Instead of using the climatological cloud fraction 
data, the range of recession coefficient due to global insolation from totally clear to totally cloudy 

t 

skies is calculated then as shown in Fig. 5 much of the scatter in the observation appear reasonable. 

One should, however, note that our modification of cloud fraction data without adjusting the 

* 

other ciimatologic parameters may not depict a realistic situation since it is known (Sellers, 1965) 
that air temperature, for example, is correlated with global insolavion and hence cloudiness. Tlie 
higher rece^on coefficient for .April-May period indicates that the rate of soil drying was slower 
than the calculated rate. During March soil remains largely snowcovered, and dormant and emerg- 
ing vegetations blanket fairly wet soil during April-May. Evaporation is expected to be slower 
and probably does not occur from entire i).3S m surface layer (assumed in this paper). Some dis- 
crepancy With observations is understandable. In the versatile water budget model (Baier and 
Robertson, 1966) the thickness of the soil layer involved in moisture dynamics is adjusted season- 
ally to take into account the changes in plant root activity and growth. 

Saxton and Lenz (1967) obser/ed that recession coefficients for the Colby watershed were 
generally higher than those for Fennimore. Since these two watersheds are only 240 km apart 
It is reasonable to assume that the climatological parameters would be about the same. Assuming 
this climatological equivalence, the difference in the recession coefficients should then arise from 
the difference in the soil texture which determines the available water, Z (9f^ - d^). From eqn. 

. (6) one would obtain 

0.7 

^Colby “ ^Fennimore (S) 
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The observed mean recession coefficient is 0.92 during July at Fennimorc, which would imply 
the con-esponding value for Colby should be 0.943. Tliis calculated value for Colby agrees well 
with the observed value of 0.95. 

4.0 CONCLUSION’. 

We discussed a formal basis for the API model and a rational method for calculating the re* 

cession coefficient. To the extent that one can use the <dimatological (or pan evaporation) data 

> 

to obtain potential evapotransphation and soil texune information, it will be possible to calculate 
a good fint approximation for the recession cy?efficient, and simulate moisture conditions in large 
areas.' Microwave remote sensing of API would be useful in updating the atmospheric general cir- 
culation models and watershed runoff forecasting. 
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APPENDIX: A PARAMETRIC MODEL FOR DAILY TOTAL INSOLATION 
Tlte insolation is cleady a m^'or factor in calculating evapotronspiration. In the absence of 
direct observations, either climatological data or a parameteric model may be used. With a para- 
' metric modei one con simulate in olation as it is affected by changes in cloud type and cloud 

cover as well as the atmospheric parameters and surface albedo, A parametric model (Choud- 

* 

hury, X9S1) for calcinating daily total global insolation is recently developed and tested against 
observations mid a set of exact radiative transfer calculations. Pertinent equations of this model 
are given below. 

Tlie daily total global insolation (KJ m"^) is 

I * 

86.4 S(,r 

S ■ [h sin ^ sin S + cos ^ cos 5 sin hj . 

IT 

T (m, c, r, a, w, u, |J, p) (A.-l) 

S(j * 1353 (W m““); solar constant 
, , 2a-(N-94) , 

r * 1 -■ 0.0335 sm [ J ; mean sun-earth distance (A-2) 

365 

N: Julian date 

h ■ 7 T - cos'* (tan ^ tan 5) (A-3) 

<f>: latitude 

5: solar declination 
1.275 

m * 0.105 + —— 7 ; — rr mean air mass (A-4) 

cos - 5) 

c: fractional cloud cover 
r: cloud optical tliickness 
Oj ; surface albedo 

w: atmospheric precipitable water (cm) 
u: ozone path length (cm) 
turbidity 

pi surface air pressure (mbar) 
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T (m, c, T, ot^, w, u, p, p) 

T(m,w, u,0,p) ( 1-c + cTgi \ 

( 0.0658 + 0 H ^ “s “cl ) 

“ “» I 0.9606 + ^ ) 

* 

; atmospheric tra/ismission function (A-S) 

. 0.97 (2 + 3/„,) 

T,i * ; cloud transmission function (A-6) 

" 4 + 0.6 f 

— ; diffpse albedo of cloud (A-7) 

^ 4 + 0.6 T 

T (m, w, u, 0, p) 

» 11 CjU (1 - Sa) (1 - Sd) + 0.5 Sa + 0.75 ] 

: clear sky atmospheric transmission function for non-reflecting 
surfaces (A-8) 

€j: absorption coefficients for water vapor (j*l), carbon dioxide 0*2), 
ozone 0*3), oxygen 0=4) and aerosols 0=5) given in Hoyt (1978) 

Sj, Sjjt scattering coefficients, respectively, for Rayleigh and aerosols given 
in Hoyt (1978) 


Gimatological data sources for the atmospheric parameters are cited in the text. Montlily 
average values of the cloud optical thickness are obtained by matching calculated and observed 
insolations. 
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FIGURE CAPTIONS 


Figure 1. Comparison of calculated amf observed (.) (dear sky insolations at Rockville, MD. Tlie 
unit for insolation is in accordance witli observations; to convert into SI units, I ly » 
41.87 KJ m’2. 

Figure 2. Calculated global insolation and net radiation for south-central Wisconsin, 

Figure 3, Comparison of simulated (.) Penman evapotranspiration and daily net radiation witli 
the regression equation from Tanner and Pelton (1960). 

Figure 4. Comparison of observed (,) and simulated annual trend of the recession coefficient 
(solid line based on Penman actuation, dashed line based on Tanner-Pelton regression 
equation). The observed data points are from Saxton and Lenz (1967). 

Figure 5. Tlie range (vertical bars) of simulated recession coefficients due purely to cloudiness 

conditions. (Penman equation). The observed data points (.) are from Saxton and Lenz 
(1967). 


TABLE 1: Climatological data for south-centnl Wisconsin Oatitude 44“N) used in the simulation of the recession coefficient. 
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Figure 1. Comparison of calculated and’observed (.) clear sky insolations at Rock- 
ville, MD. The unit for insolation is in accordance with observations; to 
convert into SI units, 1 ly = 41.87 KJ 
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Figure 2. Calculated global insolation and net radiation for south-central Wisconsin. 
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RECESSION COEFFiaENT 
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Figure 4. Comparison of observed (.) and simulated annual trend of the recession 
coefHcient (solid line based on Penman equation, dashed line based on 
Tanner-Pelton regression equation). Tlie observed data points are from 
Swton and Lenz (1967). 
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Figure 5. The range (vertical bars) of simulated recession coefficients due purely to 
cloudiness conditions (Penman equation). The observed data points (.) 
are from Saxton and Lenz (1967). 
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